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Abstract 



We investigate superconductivity in a grand canonical ensemble with fixed 
number parity (even or odd). In the low temperature limit we find small 
corrections to the BCS gap equation and dispersion E(k). The even-odd free 
energy difference in the same limit decreases linearly with temperature, in 
accordance with the behavior observed experimentally and previously arrived 
at from a quasiparticle model. The theory yields deviations from the BCS 
predictions for the specific heat, ultrasound attenuation, NSR relaxation rate, 
and electromagnetic absorption. 
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I. INTRODUCTION 



An extremely interesting recent development is the experimental study of capacitatively 
isolated superconducting metallic islands to which electrons can be added one by one. Two 
different experimental probes of such islandsB'i have demonstrated that even when the num- 
ber of electrons is extremely large — of the order of 10 9 — a periodicity is observed for every 
two electrons added. Since the microscopic BCS theory of superconductivity^ is a pairing 
theory, it is natural to assume that the observed alternation is due to the evenness or odd- 
ness of the number of electrons on the island. At first sight one might have guessed that 
adding electrons would involve a localized state near the injection site, first filling it and then 
emptying it to form a pair. However, the temperature and magnetic field dependence of the 
effect confirm that for the reported experiments the relevant energies relate to superconduct- 
ing correlations. Since these experiments probe fundamental properties of superconductors 
their precise analysis is important. What makes this problem particularly interesting from 
the point of view of theory is the fact that superconducting order is off-diagonal in number 
so that the treatment of fixed number is problematical. 

Tuominen, Hergenrother, Tighe, and TinkhamS have suggested that equilibrium proper- 
ties can be calculated from the partition functions 

Ze/o = \ III [1 + ± II [1 - j > (!•!) 

where e means even and o odd, and E^^ are the quasiparticle energies. The sign ± retains 
only configurations containing even or odd numbers of quasiparticles in a grand canonical 
ensemble. This approach has also been adopted by Lafarge et al.i in interpreting their 
experimental results. Now, Eq. (|l.l|) cannot be quite right because the BCS theory is a 
self-consistent theory in which the energy is not simply a sum of quasiparticle energies. On 
the other hand, the number parity of a state with N quasiparticles is the same as that of a 
state with N electrons. In this paper we explore the consequences of treating the islands by 
considering the appropriate generalization of (|1.1| ), namely 



-Ti{[l±(-l) N ]e-^ n -^}. (1.2) 



We obtain the free energy for both even and odd case and show that in the low temperature 
limit it reduces to the ansatz ( p. . 1|) . Using the same framework, we calculate the ultrasound 
attenuation, NSR relaxation rate, and electromagnetic absorption and find corrections to 
the standard BCS results at low temperature. 

The method of Eq. fll.lO is normally used in another context where fixed number is 
important: nuclear physicsa. There energy-levels are calculated using an unrestricted grand 
canonical ensemble but only those with the appropriate quantum number for the nucleus 
in question are kept. Since we are here interested in thermodynamic properties which are 
related to fluctuations, it would seem more correct at the outset to suppress states with the 
wrong number parity in the statistical ensemble, as in ( |1.2|) . 

In the BCS theory one fixes the average number of electrons N by adjusting the chemical 



potential. The fluctuations in particle number (of order J (NA/Ep) at T = 0, where A is 
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the energy gap and Ep the Fermi energy) are negligible for a macroscopic system. Since we 
are interested in a small free energy difference in the present work, we keep track of number 
fluctuations in both the even and the odd ensembles of ( |1.2| ), and explicitly check that they 
do not affect the physical properties we calculate. 

This paper is organized as follows: in Sec.([II|) we calculate the partition functions for 
the parity restricted ensembles; in Sec.(|T|) we minimize the grand potential and obtain 
the gap equation and dispersion for even or odd parity. We find only small corrections 
to the BCS results. Sec. ([TV]) is devoted to extracting thermodynamic properties from the 
grand potential. The differences between the results obtained with the method of Eq. ( |1 . 2| ) 
and the nuclear physics method, Eq. fll.lj) , turn out too small to be measurable with the 
experimental techniques presently available. All subsequent calculations are made within 
the approximation We evaluate the even/odd free energy difference paying special 

attention to number fluctuations. The entropy and specific heat of the even and odd parity 
systems are also calculated: the results show deviations from BCS. In Sec.(|Vp we show how 
the number parity restriction affects the response of an even- or odd-N superconductor to 
various external probes. The final Sec.flVTD contains a summary and discussion. 



II. CALCULATION OF PARTITION FUNCTION 

In principle the most suitable framework for discussing the problem would be to consider 
N particles interacting via the BCS interaction and construct a BCS-like trial ground-state 
wave function and excited states with a fixed number of particles!. For our purposes, 
however, this scheme runs into technical difficulties. The fixed-N BCS wave function cannot 
be normalized exactly, and a saddle point approximation yields the grand-canonical BCS 
results. 

If one tries! to use the even and odd ground states in the usual grand-canonical method, 
the difference between them will be practically removed by the parity-violating grand- 
canonical process of a single-particle exchange between the system and the particle reservoir. 

In the partition function (|1.2|) the projection operator [1 ± (—1)^/2] suppresses states 
with odd and even number parity, respectively. The expectation value of an operator in 
such an ensemble has the form 

(0) e/0 = Z^-^^Oe-^-^}. (2.1) 

The reduced BCS Hamiltonian in the expressions above, 

U- fiN = J2 e kc[ a c k(7 - Vfefc'4T c -^ c -fc'i c fc'T> (2.2) 

k,cr kk' 

(where the single-particle energies are measured from /i) can be transformed, using the 
Bogoliubov-Valatin transformation! 

Iko = u k c k ^ - v k cl ki , (2.3) 
Tfei = v kc{^ + UkC-ki-, (2.4) 
«* + ** = !. (2-5) 
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into the expression 



U- fiN = H Q + H I + U. (2.6) 
Here TLq corresponds to the kinetic energy of the Bogoliubov quasiparticles 

W Q = E S *(7iSo7fco + 7iti7fci)- (27) 



where E k will be obtained later on by minimizing the grand potential. Note that if 
were substituted into Q1.2| ), Eq. ( |L1| ) would result. However, there are two additional terms 
in (|2.6p . Tii describes the interaction between the elementary excitations. 



k L 



e k (u\ - vl) + 2u k v k Uk'V k/ Vkk' - E k 



(7fc 7fco + 7fci7*0 



v kk'UkV k u k 'V kl (jl j k0 + 7fci7fei)(7fc'o7fc'o + ihilk'i)- (2.8) 



k,k' 

Finally, U is a c-number term depending on temperature : 

U = 2J2 - v kk>u k v k u k >v k , . (2.9) 

h h h< 

We now make the BCS mean field approximation by assuming that the low temperature 
dynamic properties of the system are governed by the independent quasiparticles described 
by T~Cq , and that the interaction term can be replaced by its thermal average 

H - fiN ~Hq + C e / , (2.10) 

where the temperature dependent constant C e / a is defined as 

C e / Q = (Hl) e /o + U — 

+ 2^2e k v k -22V kk >u k v k u k >v k i. (2.11) 



k,k' 



In this approximation it is possible to calculate the even-odd partition function explicitly. 
The result can be given in terms of the quantities 

Z ± = Y[(l±e-^), (2.12) 

k.a 

Here f + (E k ) is the usual Fermi function. Then the even-odd partition function has the form 

Z e/o = ^(Z+±Z.)e- pc '/». (2.14) 
In the same mean field approximation C e / Q becomes 
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C e / — C + + SC e / Q , 
SC e/0 = ± Z f +± ~ z C+) - (2-15) 

The function C + in the above expression is the BCS result for the constant C, 
C + = 2 £ e k v 2 k + 2 [c*(«2 - i# - £*] /+(£*) 
- v kk>u k v k u k ,v k , [1 - f+{E k )) [1 - , (2-16) 

k,k' 

whereas C_ can be obtained from C + by a replacement of f + (E k ) with f_(E k ), 

C- = C+[f+{E k ) ^ f-{E h )]. (2.17) 
Inserting these relations into Eq. (|2.15|) we find 

xn 2 r _ [e fe (wj; - w 2) _ s fc ] ^ coth/3£ fe . 

sc */° = -T^Tt Trmz\2^ • h o^ + ^2^V kk ,u k v k u k ,v k , }■ 

1 ± T7fc, CT coth eft L Y smh P^fc smh ^ fc ' 

(2.18) 

The even-odd grand potential can be directly obtained from the even-odd partition function 
via Q e / = — ksT hi Z e / Q . Casting Q e / a in a form that makes the even-odd corrections 
explicit, we finally obtain 



^ e /o — &>bcs — ~o In 



1/r i±n tanh — 



(2.19) 



The last two terms in Eq. ( j2.19| ) represent corrections arising from even-odd effects. 
These corrections vanish for large enough temperature and/or volume. The easiest way to 
see how this happens when the temperature is increased is to calculate the product 

TJtanh^ = e £^ lntanh ^. (2.20) 

k,a ~ 



For A <C ksT, so that E k ~ je^l, one obtains 

rjtanh^ = e- 3w27V / 2 ^. (2.21) 



where N is the total number of particles and ep is the Fermi energy. Since all even-odd 
corrections are proportional to some power of this product, such effects are negligible at 
temperatures close to T c . As we will see, a much lower temperature limit for the onset of 
even-odd effects can be set if we start from the low temperature limit and calculate the 
temperature where the even-odd free energy difference vanishes. In fact, when A ^> ksT, 
the relevant excitations are E k ~ A. Defining N e ff bylll 
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rj t anh^ = ( tanh ^) ^ ( 2 - 22 ) 



AT, 



// 



k.a 

direct calculation shows that 

E 



poo hr , 

N eff ~ 4iV(0)V ^ /ff 2 _ A2 e ( " 2N(0)Vy/2nk B TA, (2.23) 

where V is the volume of the system. Since tanhx < 1 for all finite arguments x, the 
product is again negligible for large enough N e ff, i.e. large volume or particle number. For 
the samples used by Tuominen et al. and Lafarge et al. N e ff ~ 10 4 . 



III. MINIMIZATION OF THE GRAND POTENTIAL 



Before performing the minimization it is useful to rearrange the terms in Eq.( [2.19D as 
follows 

Q e/o = V, B cs(f+(E k ) -» f{Ek)e/o) ~ 4^2V w u k v k u k >v k > (n k n k ,) e/o - (n k ) e/o (n k ,) e/o , (3.1) 

k,k' 

where n k = ^\^ k is the quasiparticle number operator and 

f+(E k )Z+±f4E k )Z^ 



f(E k ) e /o = (n k ) e /o 



Z,±Z_ 



(3.2) 



Also 



n B Cs(f + (E k ) ^ f{Ek)e/o) 



Y. e kV 2 k 



k.a 



+ [ e k( U l ~ V k) ~ E k fe/o(E k ) - Vkk>UkVkU k >Vk> 1 - 2f e/o (E k ) 1 - 2/ e /o(-Bfc/)] . (3.3) 
fc,cr k,k' 



As we shall see later on, the last term in Eq. (|3.1| ) is small, but nonzero. This is an explicit 
manifestation of the fact that the thermal Wick's theorem is weakly violated in the even/odd 
ensembles. The grand potential is a functional of two independent functions, say E k and 
Vk, Q e /o = &>e/o{{E k ,v k }). Minimization with respect to these parameters will give two 
variational equations, which in the BCS case yield the energy spectrum and the gap equation. 
Following the same strategy, we obtain 



2\ A e/o 



2e k u k v k - (u k - v k )A k 
E k - e k (ul - vl) - 2u k v k A 



A 



e/o 
k i 



e/o 



B 



e/o 



Here A e J° is defined 



as 



Afc /0 = V k,k>U k >V kl [l - 2/ e / (£*/)]. 
k' 



(3.4) 
(3.5) 



(3.6) 



Note that u k and v k are also parity dependent. We have omitted the e/o labels for the sake 
of simplicity. We shall see later on that the gap in the spectrum is somewhat different from 



A 



e/o 
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Both A e k /o and B e k /o were generated by the residual, Wick's theorem violating term in 
Eq.(0- 



Af = 4 



u k~ v k 



, .... Y.V kk >u k v k u k rV kl (n k n kl ) e/o - {n k ) e / {n k ,) e/o 

1 - lj e / {E k ) k u 



(3.7) 



B 



e/o 



U k Vk 







-df e/0 (E k )/dE k 



k,k' 



y2V kk 'U k v k u k 'V k '—— {n k n k >) e / - (n k ) e / (n k >) e / 



(3.8) 



These coefficients also vanish in the thermodynamic BCS limit, since then, for reasons 
similar to those given above for the grand potential, even/odd effects are small, and Wick's 
theorem is restored. Also note that, since A e k /o and B e k /o are proportional to u\ — v\ and u k v kl 
respectively, they can be interpreted as changes in H\ k ° due to Wick's theorem violation. 



A. Exact solution of the minimization equations 

It is possible to give a formal exact solution to fl3.4| ) and ( |3.5| ) by making the usual change 
of variables u k = cos 8 k ,v k = sin 6 k . We obtain from the first equation of minimization, ( |3.4| ) 
a solution similar to the BCS case 



u k 



v k 



\|2V (el + A' e/o y/ 



\ 



+ A 2 , W 2 

1 e/o) 



but now A e / is given by 



A 



■e/o 



(n k n k/ ) e / - (n k ) e /o(n k ') e / }, 



1 — 2fe/o(E k ) 

The spectrum E k can be obtained from ( |3.5| ) by using the above results. 



x 



Ek = (et + A 2 e /o) 1/2 



+ 



A 



—2 : . Vk,k' 
A . U/2 



A 



■e/o 



(4 + A e/o ) 1/2 ^ ' 2(4 + A 2 /o )V 

9 x r 

(n k n k >) e / — (n k ) e / (n k ') e / 



■l-2f e/o (E k ) df e/o (E k )/dE k dE k 
The expressions above contain the quantity 

p e / (k,k') = (n k n k ') e / - (n k ) e / (n k >) e / . 
This can be evaluated explicitly. For k ^ k' 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



(3.13) 
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Pe/o{k, k ) = — — 

ff+(Ek)Z+ ± f-{E k )Z\ (f + {E k ,)Z + ± f-(E v )Z. 



\ Z + ±Z_ )\ Z + ±Z_ 

Inserting the expressions for f± and Z± [ Eq. (|2.13| - 12.12|) ] , we get 



(3.14) 



ys _ cosech/?E fc cosech/3E fc , 

Pe/o( ' } " (i±n felCT coth^)(i±n felCT tanh^)- (3 - 15) 

On the other hand, for k = k' p(k, k) has the usual form for noninteracting fermions 

Pe/o(k, k') = f e/o (E k ) - f e/o (E k ) 2 (3.16) 



The results obtained so far are exact. The expressions for A e k /o [ Eq. ( PIP ] and the 
dispersion E k [ Eq. (|3.12|) 1 constitute a closed system of two self-consistent equations which 
in principle can be solved iteratively. It is, however, possible to give approximate analytic 
result for the low temperature limit. 



B. Low temperature approximation of the exact solution 

The temperature range where the experimentsB! revealed even/odd effects is 50mK < 
T < T* ~ 200mK. In this temperature range, for the given island parameters N e ff ~ 10 4 . 
Then the following inequalities hold: 

exp(-/3A)< -J- < 1. (3.17) 

^eff 

Here A is the smallest value of the dispersion E k for both the even and odd cases. As we 
shall see the zero temperature limit of the gap in the even and the odd case differs, but 
by a small amount only. So in first approximation A is the same for both parities. These 
inequalities enable us to systematically keep track of the terms in an expansion of the exact 
result and to make a controlled approximation. In this regime 

1 1 



IE,, - ~" ' ( 3 " 18 ) 



1 ± Uk,a coth fifk i±(i + 2N eff e-? A «/° 
1 1 



1 ± Uk,a tanh £ffc 1 ± (1 _ 2A^ e// e- /3A "/°) 



,-rO 



(3.19) 



Here A e ^ = A e / D (T = 0) The parity dependent average occupation number f e / (E k ) then 
becomes 

, / tp \ f i tp \ cosech/3£' fc 
Je/o{E k ) = U{E k ) - — — pEk - 

1 ± flfca COth 

exp(-(3E k ) —J—. (3.20) 

1 ± (1 + 2iV e// e- /9A «/°) 
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For the even parity we have 



f e (E k ) = N eff exp[-(3(E k + A°)], (3.21) 



whereas for the odd case 



f Q (E k ) = exp(-^ fe )(l + (3.22) 
Similar calculation leads to approximate results for p e / {k 1 k'). We obtain for k ^ k' 



Pe {k, k 1 ) = exp[-(3{E k + E w )l ~ 0exp(-2/?A) (3.23) 
k>) = expl-W + ^-gg)] „ 0{VNlitY (324) 



We see that p e / Q (k, k') is smaller than the even/odd corrections coming from f e /o{E k ) by at 
least a factor of 1/N e ff, and therefore can be safely neglected in Eq. ( [3.1 1|) and Eq. ( p,12|) . 
The contribution from the case k = k! is also negligible, since it represents only one term in 
a summation over k'. Furthermore 

1 dp e / (k, k') , exph^-A^)] 



-df e/o (E k )/dE k dE k N t 



eff 



Finally, using the results we obtained above, it is possible to give a somewhat simpler, 
approximate expression for A e / Q and E k . 



Al /D = E 2 A £ , 1/2 [1 - Ve/o{E k >)] + O(-jrjS-). (3-26) 

2(4 + A e/o )V2 iV e// 



/ A 1 I 



^^(4 + a^) 1 / 2 

(4 + A^) 1 /^ 2iV e// (4 + A e >/ 2 ^/ 



Notice that within this approximation A^/° = A e ^ . From Eq. (|3.26|) we can now extrac 
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the low temperature behavior of A e / . Using E k ~ (e| + A^^) 1 / 2 , we obtain for the even 
case 

A e (T) ~ A BCS (0) [l - (4vr£; B TiV(0)V) exp(-2/3A e (0))] . (3.28) 

Here A B cs , (0) is the zero temperature limit of the BCS gap for the same interaction V k>k >. 
Note that A e (0) = Abcs(0). We used the expression ( |2.23| ) for N e ff to obtain the prefa- 
tor of the exponential. The most significant difference between the above expression and 
the usual BCS result is the exp(— 2/3 A) decrease of the gap, to be contrasted with the 
exp(— (3 A) behavior predicted by the conventional BCS theory. Furthermore the prefactor 
of the exponential is linear in T, whereas in the BCS is proportional to T 1//2 . For odd parity 
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A (T) ~ A 



BCS 



(0) 



-(^) 1/2 exp(-/3A o (0)) 



2A o (0)iV(0)V v A o (0) 
Note that in the odd case A D (T) saturates to a value that is less than the BCS gap: 



(3.29) 



A o (0) ~ A BCS (0) - 



1 



2JV(0)V 



(3.30) 



The difference, however, is very small, ~ 10 A A B cs- Finally, let us examine Eq. fl3.27p which 
gives the approximate dispersion. We can evaluate the sum in the second term and obtain 



V k , k . exp[-0(E kl - A° e/o )] 2vk B T ^ /2 
r 2( e ^ + A^ /o ) V2 - ^ Ae/O (0)' 



(3.31) 



where A = \Vk,k'\avN(0)V — 1/4 is the superconducting coupling constant. The smallest 
possible value of the dispersion, 5 e / = E(k = k F ), is slightly different from A e / . In fact, 
from Eq. (|3.27| ) we get 



5 e /o = A e / 



IT 



A , 27tk B T 



1/2 



iV e/ /A e/o (0) 



(3.32) 



This means that low temperature behavior of the spectral gap 5 e / a will be somewhat different 



from that of A 



e/o 



6 e (T) ~ A BCS {0) 



A 



2A o (0)iV(0)V 



- (47rA; i jTiV(0)V) exp(-2/3A e (0)) 



(3.33) 



S (T) ~ A BCS (0) 



1 - 



1-A 



,2ixk B T^ l/2 



2A o (0)iV(0)V v A o (0) 
Clearly, both S e (0) and S o (0) is smaller than A SC5 (0) 

A 



exp(-/?A o (0)) 



6 o (0) ~ A BCS (0) 



2N(0)V 



(3.34) 



(3.35) 



S o (0) ~ A BCS (0) 



1-A 
2iV(0)V 



(3.36) 



Nevertheless, since A < 1, the odd gap is depleted somewhat more than the even gap. The 
deviations from the BCS results found in this section are probably too small to be experi- 
mentally observable in the quantities we investigate later on. Therefore, in the subsequent 
calculations we ignore any deviation from the BCS gap equation and dispersion. 
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IV. THERMODYNAMIC PROPERTIES 



Now we can estimate the grand potential itself and the even-odd grand potential differ- 
ence in the low temperature limit k^T <C A, by using the standard BCS relations for E k 
and v k in eq. Q2.19p . In this limit 

o 

5C e / 



x 

kk 



1 ± (1 + 2N eff e-P A 



£ -=^V kk> e~^ E ^ (4.1) 



is exponentially small for both even and odd parity. In fact, after neglecting the parity 
dependence of C, the results obtained within this approximation do not differ from those 
obtained from Eq. ( |1.1|) . For the even-odd grand potential difference, calculated at a fixed 
chemical potential /i, we obtain 

Sn e/o (jJL) = Q {fi) -n e (ji) = A- k B T In N eff . (4.2) 

The experimentally relevant quantity is the free energy difference of two isolated (canonical) 
systems with 2N + 1 and 2N particles respectively. This can be related to SQ e / through 
the relation 



e 



N 



where Fpj(T) is the free energy energy for a canonical ensemble. Using a Gaussian approxi- 
mation about the peak of the summand, one obtains 

1 /AiV 2 \ 

n (ji ) - fi e (/i e ) ~ F(2N + 1) - fi (2N + 1) + F(2N) - /i e 2iV - — In (4.4) 

where AN^ Q are the number fluctuations in the grand ensembles ( |1.1| ). The investigation of 
how the number fluctuations depend on the number parity of the ensemble is an important 
test of the present model, since one must check that the fluctuations in the parity restricted 
grand canonical ensembles do not wash out the free energy difference one is attempting to 
calculate. 

The number operator written in terms of the Bogoliubov creation and annihilation op- 
erators has the form! 

N = E v l + ( u l - vl) [llolko + 7-fci7-fci] + 2u k v k [7lo7-fci - 7fco7-fei] • (4-5) 

k 

The last two terms are off diagonal in quasiparticle number, and therefore they do not 
contribute to the average particle number. ( They, however, turn out to be crucial for the 
mean square fluctuations AiV^ ). The average quasiparticle number is 

< N >e/ =X>* + (4 - vl) f+Z z +± J z Z - - (4-6) 
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Now we can fix the chemical potentials for the even/odd systems in Q4.4|) , since they satisfy 
the above equation with < N > e = 2N and < iV > = 2N + 1, respectively. From (^4.6| ) we 
obtain that p D = p e + <5/i, where <5/i oc ep/2N is of the order of the level spacing (much 
smaller than A). 

When calculating the mean square fluctuations AN 2 , Q =< N 2 > e / Q — < N >\i by using 
the operator expression in eq. (|L5|, we see that the off-diagonal terms now give a nonzero 
contribution, proportional to (A/Ek) 2 

+ 4(f + (i - /+>» - S0z + ~ (1 " f - r n- ^ 

At low temperatures, when f± — > ±exp(— (3 A), direct calculation shows that: 

^ N e/o = 4E( A /^) 2 + 0(1/N eff ) t (4.8) 

k 

regardless of the even or odd nature of the ensemble. This very important result reassures us 
that thermal and quantum fluctuations in the particle number do not dominate the difference 
between the grand canonical ensembles of even and odd parity in Eq.( [4.4|) . Therefore this 
modified grand canonical method can be used to calculate even/odd differences in canonical 
quantities. One should note that this property is intrinsic to the superconducting state: at 



T = the fluctuations are of a purely quantum nature AN e / oc ^JNA/Ep, and decrease as 
the superconducting gap A is depleted. In contrast, the same quantity for the normal state 



is AN e/o oc ^JNT/Tp -> as T -> 0. 

Expanding f2 e (/i e ) ~ ^(aO + Sfi2N on the left-hand side of Eq.( [4.4|) , and using 
and (|4.2j), we obtain : 

F(2N + 1) - F(2iV) - /i = A - k B T In N eff . (4.9) 

This expression, previously given in refs.BJl, has been used to provide a reasonable account 
of the highest temperature at which the 2e periodicity of the I-V characteristics for the 
single-electron-transistor is experimentally observed by Tuominen and collaborators^: T* = 
A/fcelniVe// ~ 300 mK. It also agrees with the temperature dependence of bT e jo measured 
by Lafarge et alM 

The method presented here allows us to calculate other possibly measurable properties 
of the superconducting state with restricted particle number parity. For example, we can 
calculate the entropy from the expression 

S e /o = -k B Tr{p e/o lnp e/o }, (4.10) 

where the density matrix is given by 

p, lo = Z : f^^e-^\ (4.11) 

After performing the parity restricted trace and using the results we obtained previously for 
Z ejo ( Eq. (|JJD ) we get 
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S e /o = S B cs + 7f,^ E kbfe/o + k B hi | (l i^anli^" 1 



T 

where Sbcs is the entropy of a BCS superconductor 



2 

fe,. 



(4.12) 



Sbcs = -k B Eft 1 " /+) ln (l " /+) + /+ ln f+h (^) 

k.cr 



and 5/ e / is the even/odd change in the quasiparticle occupation /. 

U-f- i i 

1 ± Z + /Z- sinh pE k 1 ± n fcCT coth £ffc 

In the low temperature limit the entropy for the even-parity system is 

S e ~ k B (3AN 2 eff e- 2f3A , (4.15) 

whereas for the odd-parity system we have 

S°~k B hiN eff + S BC s. (4.16) 

Here Sbcs = kBN e ff((3A)e~ f3A . The specific heatEl can be obtained directly from the above 
expression via 

ace/o 

C e/o = T _^_ (417) 

More specifically, for the even case 

C e v ~k B (N eff (3A) 2 e- 2 P A . (4.18) 

Note the exp(— 2/3 A) dependence of the heat capacity for this case, different from the usual 
exp(— (3 A) in the BCS case. This discrepancy is the direct consequence of the number parity 
restriction applied to the system. The odd-parity system has a specific heat that is finite at 
all temperatures due to the presence of an unpaired electron: 

C'-IfcB + Cf 75 . (4.19) 

Finally, we would like to make some technical remarks regarding the calculations pre- 
sented in this paragraph. The result in Eq. ( |4.6j ) can be obtained in another way as well, by 
using the thermodynamic relation < N > e / = —dQ e /o/dfi keeping in mind the requirements 
from minimization, dQ e /o/dEk = and dfl e /o/dvk = . In contrast to this, the thermody- 
namic relationship AN 2 j = jjd < N > /dfi does not give the correct result for the mean 
square fluctuations: the crucial off-diagonal terms do not contribute to < N > and as a re- 
sult are missed in the expression for AN 2 ^ derived in this way. This reflects the peculiarity 
of the superconducting state. Yet another possible source of error is encountered if one tries 
to calculate the entropy by directly applying the thermodynamic relation S = —(dQ/dT)v 
instead of using the expression (|4.10|) . The reason is the following: the temperature depen- 
dence of the function C =< Tlj > is an artifact of the mean field approximation. Therefore, 
when calculating the entropy in this way, one should not take the derivative of C. In fact, 
by evaluating S e / with the help of Eq.( f4.10| ), one can explicitly show that within the mean 



field approximation all contributions from C e / Q drop out. Exactly the same situation arises 
in the ordinary BCS calculation. 
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V. TRANSPORT PROPERTIES 



Here we calculate the response of the parity restricted superconducting state to various 
external probes. The acoustic attenuation rate of ultrasound in a superconductor, within 
the framework of the BCS theory!, is given by the expression 

= goggO r dE / E E ' ( ! _ [f+(E) _ f+(E% (5.1) 

V ; hjj Ja y/E 2 -A 2 J( E ') 2 - A 2 V EE'J lJ + K ' J + K >h 1 ; 



where E' = E + Two. In the parity restricted ensemble the thermal factor is replaced by 

/tW _ /t(gH MlAg |^g)^ (52) 

As a result, we obtain for the ratio of the rates in the superconducting and the normal phase 

(^(T))e/ = 2/(A) _ l±n fei(T coth^sinh/5A- (5 ' 3) 

In the low temperature limit, for the even case the ratio drops faster than in the BCS case 

(as(T) 



-ol n {T) 



) =2N eff e-^ A . (5.4) 



This should be contrasted with the exp(— f3A) dependence of the BCS ratio. The odd case 
also shows a deviation 

^a N {T)So N ef f 

Both terms are of the same order of magnitude at T*, the onset temperature of even-odd 
effects. As we will see later on the low temperature response functions to other exter- 
nal probes, like the NSR relaxation rate and electromagnetic absorption, will have similar 
behavior as the ultrasound attenuation calculated above. 
The BCS result for the nuclear spin relaxation rate is! 

R S (T) _ r [EE' + A 2 ]/ + (£)(l - U(E') 

R N (T) Ja [E 2 - A 2 ]V2[(£')2 _ A 2 ]V2 " 1 ' > 

Within the even/odd ensemble, the thermal factor is modified to 
f , W1 f f+(E)Z + ±f-(E)Z- U(E)U(E')Z + ±f4E)f4E')Z_ 

f + (E)(i-f + (E))^ z^TzZ z^TzZ • ( } 

An explicit estimation shows that only the first term will contributes significantly, and the 
second can be neglected. Denoting with F e / Q (E) the low temperature, E — > E' limit of the 
new parity dependent thermal factor, we obtain 

F e (E) = N eff e-^ E+A \ (5.8) 
F (E) = e-? E + e —- . (5.9) 
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Using this function, the ratio of the relaxation rates can be written as 



/ Rq(T) n r°° \E 2 + A 2 le- /3(i? - A) 

where the integral now is only a weak function of temperature. All the important temper- 
ature and parity dependence is embodied in F e / Q (A). Also note that F e / (A) gives the low 
temperature behavior of the acoustic attenuation as well, since (asf&^e/o — 2F e / (A). 

Finally, let us discuss the infrared absorption. At very low temperatures there are only 
a few thermally excited quasiparticles available, and therefore substantial contribution to 
the absorption rate is given by the pair breaking process in which an incoming photon 
breaks a Cooper pair and creates two quasiparticles. In the case of the odd-parity system, 
however, there is always an unpaired quasiparticle present, which will lead to nonzero subgap 
absorption. Let us therefore consider the case when hu < 2 A, so that only single particle 
processes are allowed. We examine the real part of the conductivity - proportional to the 
rate of absorption - which has the following BCS formi'0 

-»r /+(A ) + A' r V+{ T-£ l \ (5,1, 



a 1N (T) L ' J a £ 2 -A 2 
A calculation similar to that presented above for the nuclear spin relaxation gives 

The common feature of all these results obtained for the response of the parity-restricted 
superconductor to various probes is, that for low temperatures and energy transfers the 
exp(— /3A) dependence - characteristic of a BCS superconductor - is replaced by a more 
rapidly decaying exp(— 2/3 A) behavior for the even case, and the appearance of a very 
small (0(1/N e ff)), but nonzero background due to a single, unpaired quasiparticle, that 
can scatter, absorb, spin relax etc. Pictorially, the exp(— 2/3 A) dependence comes from 
the exclusion of the parity-violating grand-canonical processes that involve a single particle 
exchange between the system and a particle bath (excitation energy A), as well as from the 
important contribution of the momentum-conserving excitation process of a pair (excitation 
energy 2A) , in contrast to pair breaking with no momentum conservation (excitation energy 
A per particle). 



VI. CONCLUSIONS 

In conclusion, we have presented a theoretical method to investigate number parity effects 
in the superconducting state by using a grand canonical ensemble in which the parity of the 
system is even or odd. This corresponds to a permeable wall between the superconductor 
and the particle bath that allows particle exchange only in pairs. We obtained the partition 
function which in turn yielded the grand potential, entropy and specific heat of the system. 
While the variational equations from the minimization of the grand potential gave essentially 
BCS results, the free energy difference between even and odd parity states obtained is 
in accordance with the experimental and theoretical findings of Tuominen et al. and the 
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measurements of Lafarge et al.Era. The specific heat shows a more rapid exp(— 2/3A) decay at 
low temperatures; the same quantity in the odd-parity system has a finite zero temperature 
limit ks/2, in contrast to the usual BCS result. 

A simple and systematic low temperature approximation was made possible by the fact 
that for the temperature range and size of the islands investigated so far0i N e ff ^> 1 and 
therefore an expansion was possible in powers of l/N e ff. In principle, this expansion breaks 
down at extremely low temperatures, since 1/N e ff oc T -1 / 2 . In order to observe such effects 
in the temperature range presently investigated (T > 50 inK), the volume of the island 
must be four orders of magnitude less, about 10~ 19 cm 3 , which might be very difficult to 
achieve experimentally. 

We have not addressed the problem of the magnetic field dependence of the even-odd 
free energy difference. The measurements of Tuominen et alx! and Esteve et al. seem to 
suggest that the experimental results can be explained if one replaces in eq. Q4.2|) the gap 
in zero field A with the gap edge A(H) of the field dependent density of states provided 
by the Abrikosov - Gorkov theoryEl. We believe that this replacement, although possibly 
correct, is not justified well enough. The present method cannot be easily generalized to the 
Green's function technique, due to the lack of a Wick's theorem: for the parity-conserving 
average defined in eq. ( |2.1| ), (n^n^) 7^ (hk)(nk>) for a system of noninteracting fermions, 
where n& is the number operator (cf. Eq. (|3.15|) ). There is, however, a Wick's theorem for 
the two channels, with projection operator 1/2 and ±(—1)^/2, respectively, which add up 
to the parity-restricted average. Because of this, the normal-state properties and the Cooper 
instability need further investigation, to be performed in the future. 
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